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We study theoretically the responsivity of optical modulators. For the case of linear response we 
find using perturbation theory an upper bound imposed upon the responsivity. For the case of two 
mode modulator we find a lower bound imposed upon the optical path required for achieving full 
modulation when the maximum birefringence strength is given. 
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Introduction - Optical modulators are devices of great 
importance for telecom and other fields. These devices 
allow controlling the transmission T (0 < T < 1) be- 
tween input and output ports by applying some external 
perturbation. One of the key characterization of optical 
modulators is their responsivity, namely the dependence 
of T on the applied external perturbation. Enhancing 
the responsivity is highly desirable for many applications. 
This raises the question what is the largest possible re- 
sponsivity that can be achieved for a given perturbation 
mechanism employed. Here we consider this question for 
the case of linear modulators. We show that the linearity 
of such devices imposes upper bound on the responsivity. 

Perturbation Theory - Consider an optical modulator 
consisting of an optical path of length As = si — sq. 
Here we consider the case where the light passes the op- 
tical path only once (contrary to the case of a resonator 
where multiple reflections occur). At each point s along 
the optical path the field is expanded using some local 
orthonormal basis. Using the Dirac notation {bra and 
ket) U the field at point s is denoted as |V' (s)) , which 
represents a column vector of amplitudes. The equation 
of motion is given by 

^IV') = ^/c|V'), (1) 

where the Hermitian linear operator K, is the Hamiltonian 
of the system. Consider the effect of adding a small per- 
turbation eJCi (s) to the unperturbed Hamiltonian JCq, 
namely 

/C (s) = /Co (s) +e/Ci (s) , (2) 

where |e| << 1 is a small real parameter. For any given 
e the final state |V'/) = IV' (si)) is related to the initial 
state IV'i) = \tp (so)) by the relation 

|V/) = C/(e)|V.), (3) 

where U (e) is the s evolution operator for the Hamilto- 
nian JC = /Co+e/Ci from s = sq to s = si. The final state 



V'/) is filtered by a polarizer having a normalized state 
\ipp). The transmission of the modulator is given by 

T{e) = \{i,,\i;j{e))\\ (4) 

Given the perturbed and unperturbed final states, 
\ipf (e)) and \tpf (0)) respectively, what is the optimum 
choice of a normalized {ipp) that will maximize \dT/de\ ? 
Define the density operator 

p(£) = |V'/(£))(^/(£)|, (5) 
and the operator 

Ap = p{e)^p{0). (6) 
For a small e one has 



— ^-(^plAplV'j,). (7) 

The operator Ap is Hermitian, thus \tpp) that will max- 
imize \dT/dX\ is an eigenvector of Ap with the largest 
eigenvalue in absolute value. The eigenvalues of Ap are 

given by ± yT^KvTfe^KTwyF. thus 

^|<^\/l- 1(^.1 C/t iO)Uie)\^.)f. (8) 

Using perturbation expansion Q one finds to 2nd order 
in e 

{i:,\uH0)U{e)\^/M) (9) 
= 1 + iejds' (/Ci (s')) - e^Jds'jds" (/Ci (s') /Ci (s")) , 

So So So 

where the symbol () represents expectation value, 
{A{s)) — {^pilAff Itpi) for a general operator A, where 
Ah is defined as 
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Ah (s) = ul {s, So) Auo {s, sq) 



(10) 
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and Mo {s, Sq) is the s evolution operator from sq to s 
generated by ICq. Since ICi (s) is Hermitian one finds to 
lowest order in e 



= 1 - e^Jds'Jds" (A/Ci is') A/Ci (s")) , 

So So 

where A/Ci (s) = /Ci (s) - (/Ci (s)). Thus 



Jds'Jds" (A/Ci (s') A/Ci (s")) 

So So 



(11) 



dT 


1 




< 


de 





(12) 



This upper bound imposed upon the responsivity can 
be further simplified by employing the Schwartz inequal- 
ity 



is bounded by \eKi (s')| < k„ 
only for 



As > 



full modulation occurs 



(18) 



Examples - As a simple example, consider a modulator 
based on an optical fiber. Circularly polarized light is 
injected into the fiber and a polarizer located at the fiber 
end allows transmission of only linearly polarized light. 
Modulation is achieved by applying linear birefringence 
along some section of the fiber having length As. 

For the present example we chose It/;) = 1+; 2) at s = 
(|±; u), with u being a unit vector, denotes an eigenvec- 
tor of (T-u with an eigenvalue ±1), and the polarizer state 
is \ij}p) = |+; 3). Moreover, /Cq = and Ki = n-a, where 
K = (1/2) (fci, 0, 0). Integrating the equation of motion 
yields 



dT_ 

de 



<Jds\ ([AlC, is')Y 



(13) 



Two-mode Case - Consider the case where the dimen- 
sionality of lip (s)) is two. Ignoring the common phase 
factor, the Hermitian operator ICi can be assumed to be 
traceless, thus it can be expressed as 



(14) 



where Ki= \ki \ ki is a three-dimensional real vector with 
length |ki| (ki is a unit vector) and the components of 
the Pauli matrix vector a are given by: 



1 

1 



,<^2 = 



-i 

1 



,CT3 = 



1 

-1 



(15) 



It is straightforward to show that Eq. E| for the 
present case yields 



dT 

de 



<Jds'\K^ (s')| 



(16) 



Similar upper bound can be found for the an- 
gle 6 between the polarization unit vectors p (e) — 
(V-d C/t (e) aU ie) and p (0) ^ C/t (0) aU (0) 
on the Bloch sphere. Using ^2 and assuming the case 
6 << 1 one finds 



e < 2 fds' \eKi (s') 



(17) 



Full modulation between T = and T = I requires 
that the total change in 9 exceeds tt (assuming \ipp) is 
kept fixed). Thus, if the applied birefringence strength 



T (e) = sin^ 



efcisi 



(19) 



Thus at £ = the derivative \dT/de\ approaches the 
upper bound given in 1161 Moreover, full modulation is 
obtained for sq = — 7r/2efci, and si = 7r/2efci, thus also 
for this case the upper bound given by El is achieved. 

The next example deals with a modulator based on a 
transition between adiabatic to non-adiabatic regimes, as 
in Ref. 0. Consider the case where K. = k - a, 



K{s)^^(^0,^y-{jsf 



"fS 



(20) 



where 7 is a real constant with dimensionality of 
1 /length, A is a non-negative dimensionless real parame- 
ter, and |7s| < A. 

Consider the case where for sq = —A/7 the state of 
the system is a local eigenstate of /C (s) with positive 
eigenvalue, namely \ip {sq)) = | — ;3). When A >> 1 
the state evolves adiabatically Q and remains parallel to 
K (s). The polarizer is located at si = A/7 and its state 
is given by \tjjp) = l"!^). Thus in the adiabatic limit 
T — 0. Approximation solution for the case A > 1 can be 
found by considering the lowest order correction to the 
adiabatic limit , Q • The value of T is the probability 
of Zener transition to occur which can be calculated to 
lowest order 



T 



(for A > 1) 



(21) 



This approximation is compared with the calculated 
value of T (A) obtained from numerical integration of the 
equation of motion. The case A = 5 is presented in 
Fig. ^as an example. Figure[21(a) shows the calculated 
value of T (A) in the range < A < 5. Comparing 
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Comparison between the numerically calculated 
\dT/d\\ and the above upper bound is seen in Fig. |21 
(b). Contrary to the previous example, in this case the 
upper bound is not reached for any value of A. How- 
ever, in the transition region, between the adiabatic and 
non-adiabatic limits, near A = 0.695 the responsivity is 
only some 2% below the upper bound. Similarly, for the 
modulator discussed in Ref. ^^j , it was found that largest 
responsivity is obtained in the transition region between 
adiabatic and non-adiabatic limits. 



FIG. 1; Example of numerical integration of the equation of 
motion for the case A = 5. On the left the curve k, (s) is shown 
and on the right the evolution of the polarization vector p (s) 
on the Bloch sphere is depicted. 



the approximated result in Eq. 1211 with the numerical 
solution shows, as expected, good agreement for A > 1. 

On the other hand, the upper bound given by Eq. 1161 
for this case reads 



Note that the bounds discussed in this paper can be 
employed for other linear systems. For example, the 
same analysis may lead to a lower bound imposed upon 
the time required for performing a given quantum gate 
on a system of quantum bits in a quantum computer, 
when the maximum perturbation strength is given. 



dT 



dX 



< 



->^o/l 



7^0 



= TrAn 
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FIG. 2: Calculated and upper bound of responsivity. (a) 
Numerical calculation of T vs. A. (b) Comparison between 
the calculater \dT/dX\ and the upper bound given by Eq. 1221 



